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SUrtMRV 


The plenc motion of an articulated pipe made of tMO segments Is examined 
and the flow velocity at which flutter manifests Itself Is sought. The pres- 
sure In the reservoir feeding the pipe Is kept constant. In contrast to pre- 
vious works, the flow velocity Is not taken as a pre*'.rlbed parameter of the 
system but Is left to follow the laws of motion. This approach requires a 
nonlinear formulation of the problem and the equations of motion are solved 
using Krylov-Bogol luboy's method. A graph of the amplitude of the limit 
cycles, as a function of the fluid-system mass ratio. Is presented and con- 
clusions are drawn as to the necessity of considering noni {near i ties in the 
analysis. 



INTRODUCTION 


A 

A rttvIeM of existing knowledge [1j on the behavior of pipes conveying 
fluids reveals tKet all earlier analyses assume' that the fluid velocity rela~ 
tive nci tl.* ^'*pe Is a known quantity and Is unaffected by the motion of the 
pipe. In- tNese ^rev4ous works, the flow velocity Is assumed cither constant 
or period z with n prescribed amplitude and frequency. This approach elimi- 
nates the n^d to find the flow aquations of motion, Is adequate for Infini- 
tesimal trans'^'er^'' amplitudes of motion of tne pipe system, but Is Incapable- 
of predicting what will be the effect of larger amplutudes. This last short- 
coming may be of Importance when flow velocities are near critical velocities, 
that Is, velocities at which che system begins to flutter. 

It is therefore the purpose of the present study to Investigate In greater 
detail the dynamic behavior of pipes In the vicinity of critical velocities. 
Such an analysis requires a nonlinear formulation of the problem. It will be 
shown that nonlinear terms are generated by the axial motion of the pipe and 
by the flow velocity fluctuation. Thus, the relative velocity of the fluid 
will no longer be assumed a given parameter in the problem, but will have to 
be determined from its own equations of motion which will be coupled, through 
the nonlinear terms, to the transverse equation of motion of the pipe, "he 
specific system (see Fig. I) considered here is a vertically hanging articu- 
lated pipe made of two segments. Each articulation Is assumed perfect in the 
sense that it has no damping or restoring force associated with it. The fluid 
entering the pipe comes from a reservoir maintained at constant pressure and 
after Its passage through the pipe, it is discharged, tangentially to the end 
of the pipe, to the atmosphere. This articulated pipe system was selected, 


* Numbers in brackets refer to referances at the end of the paper. 


*ather th«n a continuous pipe, in order to obtain an analytical solution and 
therefore 9 et ntore complete insight into the probl€>in. The related problems of 
an articulated tv»o-iink pipe fed by a reservoir in which the pressure fluctu- 
ates with time and the problem of a continuous cantilevered pipe fid by •. con- 
stant pressure reservoir are also being analyzed for presentation in subse- 
quent papers. 

The basic assumptions made are: I) The fluid is incompressible; this is 

a reasonable approximation for a liquid since the characteristic time taken 
by a wave to travel the length of the pipe Is very much shorter than the per- 
iod of the transverse motion of the pipe; 2) The diameter of the pipe is 
small compared to its length so that the rotatory inertia of an element of 
the pipe is negligible; 3) The velocity profile of the fluid, at a given 
cross-section is uniform, modeling fairly well a completely turbulent flew; 
k) The motion of the articulated pipes takes place in a plane. 

EQUATIONS OF MOTION 

The derivation of the relevant equations of motion can readily be ob- 
tained by the direct application of Newton's law to our system (see Fig. 1). 

The first equation is obtained as a moment condition with respect to the upper 
pivot point on the free-body diagram of the two segments of the pipe. Sim- 
ilarly, the second equation Is obtained as a moment condition with respect to 
the pivot point on the free-body diagram of the lower segment. Finally the 
equation governing the motion of the fluid in the pipe is obtained by the force 
condition in the tangential direction on a fluid element and by a subsequent 
integration over the length of the pipe. The three equations thus obtained are 
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'1 \ l^l^2*2 

(ttij + '*’ ”’2^1 ” i”” ®2 
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•*■ nijij£2'' -‘^^02 ” '*' ”’ 1 '^ AjSln(02 - 0^) ■ 0 


( 1 ) 


A ® 2 

( 1 / 2 ) (m, + m2)gl2S»n02 + 0/3) (m^ + m2)A2^^2 * "’l^®2^2 


(2) 

0 2*2 
+ (1/2)(m, + m2)A^£2^0,cos(02 ' 9,) + 0/2)(m, + 0, sfn(02 - 0,) - 0 

O *• 

- Lv + tn|g (i-^cos0^ + £ 2 ^ 0502 ) “ ni^v(£j + ^ 2 ) 

* 2 2 

+ m|0| [l^ /2 + £^£ 2 ^ 05(02 ■ 0|)1 “ m^0j£j£2Sln(02 “ 0j) 

+ (i/ 2 )mj 02 £ 2 ^ - 0 ( 3 ) 


where and m 2 are the mass, per unit length, of the fluid and the pipe, 

respectively, g is the acceleration due to gravity, £^ and £2 are the 

lengths of the upper and lower segments, v is the flow velocity, is the 

. 2 

force, due to pressure, acting on the fluid at x ■ 0 and Lv represents 
the friction force between the fluid and the pipe. Fluid dynamic studies [2] 
have shown that L depends on the size and roughness of the pipe, on the den- 
sity of the liquid ar.d on I^eynolds number if the latter is not too high. For 
fully turbulent flow (high Reynolds number), L is independent of Reynolds 


3 


number and consequently is independent of viscosity. 

* 

We need noM to relate the pressure p in the reservoir to the force 
"p^" . From Bernoulli's equation, applied between the surface of the water in 
the reservoir and its outlet, we obtain: 





where the subscript "o" refers to quantities measured at the entrance of the 
pipe, h is the height of the fluid in the reservoir, p is the density of 
the fluid and p* is the pressure in the reservoir. Eq. (i») is now combined 
with Eq. (3) and after the introduction o^ the foliowing non-dimcntional 
variables; 



Eqs. (1,2,3) becooie: 


(a^ + a^/3)e^ + (a/2)^-COS (02 ~ 9j) Yd^U0| + 2aU92COs{02 " 9|) 
{a * a^/2)sin0^ YaU^si'’(02 ~ 

+ Y'jii sin(02 " 0j) " (1/2)02 " ^1^ " ^ 

k 


(5) 



O/2)s;n02 + (1/3)02 ^ l/2)a[9 jCos (0^^ - 0^ + 0,^sin(02 - 0j)j - 0 (6) 

P* - 0^[(l/2) + L^l + H* + a cos 0^ + cos 0^ - (a + l)U 
- -[(l/2)a^0,^ + o0j^cos(02 “ 6^) - a0^sln(02 ’ 0,) + (1/2)02^1 (7) 


Since we do not want here to solve the complete nonlinear set of equations, 
but only to determine the effect of the leading nonl inearltles on the system, 
we will omit nonlinear terms of order 4 ar.o higher. Also, to facilitate 
manipulations by making certain terms linear In U , we will assume that: 


U ■ LI + AU where U 
o o 


a 1 


1/2 + L 


is the steady»st?‘ 'low velocity under pressure p when the pipe is 

vertically at rest»and AU is a small fluctuation of the flow velocity such 

that 0^ - U ^ + 2U AU . 

o o 

Under these modifications, Eqs. (5,6,7) become: 


a.^ cVb 
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ya^U 

o 

2YaU 
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1/3. 
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o 
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0 

1/2 
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. ''2 


( 8 ) 
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where 


Fj - 02u(02 - e,)^/4 ©2^0(02 ' 3,)/2 - 0,ya^AU - h^2yc>tU 

+ + (o + a^/2)0j^/6 - yo>AU{02 “ 9|) 

+ YaU^^(e, - 0,)V6 - 2 yoU ALl(0, - 0.) 

O 4 I 0/1 

r . - 2^/12 - yAU 02 + 00^(02 - 0 ,)^/^ - a 0,'(02 - 0,)/2 

(a + OaU + (1 + 2L )U AU - -ae,^/2 - 0,^/2 + 0,^(a •♦• a^/2) 6,^/2 

o o I / I / 

- a0^(92 - e,) (9) 

The linear part of Eq. (C) consistsof two coupled equations with constant co- 
efficients The sciffness and damping matrix are nonsymmctric, indicating a 
nonconservative system. Eq. (9) is uncoupled from Eq. (8) in its linear pa"t 
but not in its nonlinear t-.rms. 

METHOD OF SOLUTION 

To solve the above system of equations we ./ill proceed as follows: a) 

Solve the linear part of Eq. (8) and determine 9|(t) , 02(t) for 
critical; b) Use the solution found in a) to calculate the nonlinear terms of 
(9); c) Solve Eq. (9) treating the nonlinear terms as forcing functions 
and obtain AU(t); d) Use AU(t) in the right-hand side of Eq. (8) and apply 
Kry lov-Bogol i ubov' s averaging technique. 

To find the solution to the linear part of Eq. (8) we assume a solution 
of the form: 


6 





•fid fvwttltutc tnto Eqi. (8) which yields 


(a/2)X^ 

(o* + a^/i)X^ + Yo^U X 

o 

♦ (o + - YciU ^ 

o 



- 0 


(10) 


For a nontrivial solution to exist, the determinant of the ohove matrix has to 
be zero. This leads to a lith order polynon'ial in \ . Since we (^re only 
interested in critical values of a , we know that it has to be a pure imagin- 
ary number and since the coefficients of the polynomial are all real, we can 
separate the polynomial Into Its real and Imaginary pares In the following manner 


X(X^U ayd a)/3 + U y( 1 a) - U V] - 0 
o o o 


(11) 


X^o(a ♦ 3A)/9 + X^[(o + 1)^/3 - o^'6 - uA(a/2 ♦ 1/3 - ya)] 

o 

♦ (1 + a / 2)/2 - yaU ^/2 - 0 ( 12 ) 

o 

2 

Tne elimination J X between these two equations yields: 

• I| *2 

AU-BU+C-0 
c c 


where 
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03 ) 


U i U critical 
c o 


-B ± '/b^ - kkC 

n 


and 


A 


Y^(tt ♦ 3M) . 3Y^(g/2 1/3 - Ya) 

a(I a)^ o(1 + o) 


B -- 


2Y(g 3 A) . 3Y(g/2 ♦ 1/3 ■ Yg) _ Y(g D , y , gY 


g(l a) 


2 2(1 - g) 


» Vk . (g 1)^ , 0 g/2) , g 

a g 5 2 


The substitution of Into Eq. (II) permits us to find the two complex 

conjugate, critical, eigenvalues of our system. Eq. (13) gives us two values 
for . The value of Interest Is the higher one since the other value cor- 
responds to the point where the system regains stability (see Fig. 2). To 

a 

obtain the two corresponding eigenvectors, we substitute and the eigen- 

value into either one of Eqs. (10), say the first one: 


^2j 


X.V3 + X yU, + >/2 

- J 

Xj g/2 


j - 


Since X^ and X 2 are complex conjugates, and are also complex con- 

jugates. We will not Investigate the other two eigenvalues since a previous 
analysis bv Bohn and Herrmann (3l shorirs that tney possess a negative real part and 
consequently they contribute but little to the solution. A typical root locus 
(as a function of ) for the system Is shown In Fig. 2. The motion of the 
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pipe In Its fijtter mode Is therefoi 


8^ ■ RelR|62|(cos ut + I sin wt) ♦ ^^022 (cos ut - I sin ut)] ( 1 ^) 


«2 ■ Re[02^(cos wt + I sin wt) + 022(cos wt - I sin ut)] (15) 

where lu ■ X . These equations can be expressed more conveniently as: 

0 . ■ 0 '^a^ + sln(tut + ♦.+♦) where (a + Ib) ■ R, ( 16 ) 

I 10 I 


02 ■ 0 sln(ut + ^^) 


whore 9 


A 


Re0 


21 


Rc022 ) 


+■ (“Im 02 ^ + 


Im022) 


tg^ 


o 



tg^, 


(Re02^ + Re022) 

(-!m0^^ + Jm©22) 


( 17 ) 


0 and are determined by the initial conditions. We note, in the above, 

that 021 and 022 are not completely determined since we have two equations 
and four unknown. If we want, but there is no obligation to do so, we can re- 
quire the Imaginary part of Eqs. (Ik) and ( 15 ) to be zero, which leads to: 

®21 • ®22 • 

A graphical representati .n of th>B mode shape through a cycle Is presented 
in Fig. 3 . The angular velocities of the two segments are indicated by curved 
arrows and each of the four sketches is taken at the instant when one segment 
has no angular velocity. sketches are separated by almost exactly 9P*. 
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Kow th«t wft have datamtnad th« mod# shape of Interest, we will eliminate the 
danped sisde, thus reducing our fourth-order system to second order. To 
achieve this ellalnatlon we follow Foss' work I^i] on the uncoupling of the 
equations of sot Ion of a damped linear system. His approach must, however, be 
modiflea since he assumes that all the matrices of the system are symmetric. 

To take care of this difficulty we must have recourse co the concept of the 
adjoint problem, that Is, the problem obtained by transposing all the matrices 
of our system. It can be shown that the eigenvalues of this new problem are 
the same as those of the original problem but that the eigenvectors are dif- 
ferent. Yet, these eigenvectors have the crucial property of being orthogonal 
to the eigenvectors of the original system, thus permitting the separation of 
the modes of our system. The details of how to proceed will not be given here 
because of limitations in space; they can be found In reference [S] • 

Thus, following the method described above, Eq. (8) can be rewritten as: 

i* * B^)x - 2Bi + * ‘ -a + "{ ■ * (- f u - 8) + f 7 (18) 


where 


X ■ 0/2 sln\|) , + 


(19) 


(B + 1(d) - X 


(r + Is) ■ 2 (dI 


J 


+ a^/3 a/l"\ 


a/2 


1/3, 


+ yU 


f. 7n2 


Lo 



( 20 ) 
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f • F| and F^ are dafinad In Eq. ( 8 ) 

9 ■ 8Fj 

0/3)X^ ♦ + 1/2 

(a ♦ lb) ■ F ■ - y ( 21 ) 

X a/2 

X^(a^ + a^/3) + X(Ya*0^) ♦ 3a/2 - yM ^ 

(c •*• Id) - R X 2. ( 22 ) 

• X a/2 

The damping coefficient B Is a small quantity, since wo are Interested only 
In flow v«llcltles close to critical and therefore It will be eliminated In 
Eq. (I 8 ) except In the linear damping term. Also, because of the smallness of 
that term. It has been transferred to the right-hand side of the equation, 
which contains the small nonlinear terms. The nonlinear terms F^ and F^ 
are not completely known yet because they depend on AU and A(7 . To deter- 
mine these fluctuations of the flow velocity we combine Fqs. (9) and (16) to 
obtain 

AU(a + I) + aO(I + 2L^) - (eV2)[(A + 9) - (A - B)cos 2i(» + C sin 2^] ( 23 ) 

where A ■ -a^a/2 - 1/2 + (o + a/2)u)^b^ + aa)^(a - b^) 

6 ■ -ob ^/2 + (a •*> a^/ 2 )(D^a^ + u ^/2 - ow^b^ 

2 2 2 

C ■ -oab - (2a + a )ui ab aw (b - 2ab) 

The steady-state homogeneous solution is AU •• 0 . The particular solution Is 
assumed to be of the form 


II 



(2M 


AU - 0^(K + U|Sln * UjCOS 2\{/) 

and is substituted Into Eq. (23). After a contparlson of the coefficients of 
similar terms we obtain: 


A B 

2U^(1 2L ) 
o o 


(25) 


U^(1 + 2LJC/2 " (A - B)u)(a ♦ 1) 

0 ^(1 + 2L )^+ Au^(a + 1)^ 

o o 


( 26 ) 


U 


2 


-U^O + 2L^)(A - B)/2 - Co)(a + 1) 

U ^0 + 2L )^ + 4o)^(a + 1)^ 
o o 


(27) 


Now that the right-hand side of Eq. (18) is known, we will solve this equation 
using the Krylov-Bogol iubov (K.B) method [6], keeping only the first term in 
the asymptotic expression. This method Is also known as the method of 
"averaging". 

The K.B. method essentially assumes that in the region close to the crit- 
ical point, the motion will be oscillatory, as it is at the critical point, 
except that the amplitude and phase of the motion will change slowly due to 
the effect of the small damping and of the small nonl inear i ties. 

In our particular case, the oscillatory motion of Interest is: 


X - (0/2)sin ^ - (0/2)sIn(wt ♦) (28) 

X ■ (0/2)u cos ij) ■ (0/2)w cos(u)t + ♦) (29) 
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which sctisfles che left-hand side of Eq. (18). We now assume that the small 
terms of the right-hand side of Eq. (I8) will cause 9 and ^ to vary slowly 
with time. For this to be true, Eq. (29) requires that: 

(e/2)sln t ♦ ie/2)^ cos ♦ - 0 (30) 

We also obtain, by a differentiation^ of E^.(2$): 

X • (6/2)u> cos ♦ - (0/2)u(u ♦ ^)$ln i|> (31) 

Equations (28) and (31) transform equation (l8) Into: 

(0/2)u) cos - (6/2) #u) sin ^ ■ 26x ( *~ j ^ )“ ^ ” (32) 

' r ♦ s^' s 'r^ ♦ s^ s 

which can be combined with Eq. (30) to obtain: 


a 

0 







(33) 


©♦ - 


•f 

r ♦ s' 


sr+s' s 


sin tj; 


(3M) 


Since 0 and 4 change very slowly, we assume them to be constant for tha 
duration of a cycle and integrate Eqs. (33) and (3^) to find the average of 
0 and ^ : 


e 

3 

av 


2ir o 


2ir/w 


2Bi ^ (4^) * 

r" i 3 s r ♦ s 


fu 

s 


cos dt 


(35) 
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sin f dt 


(36) 


< ■ - (4^) * f 

2iro «[ r +s's r*^s 2 


Th« danping frm, 20x , In tha argumant of tha abova tntagrals depends on 0 
to the first poMor, but tanr.i In f and g depend on 6 to the third poMer, 
which win permit us to find equilibrium points for certain amplitudes. If 
we are exactly at the critical point* B ■ 0 and, therefore, the nonlinear 
terms will decide whether or not the criticsl point Is stable. 

The evaluation of the right-hand side of Eqs. (35) end (36) turns out to 

be extremely tedious since r , s , f and g are complicated expressions. 

The details of this evaluation, which are given In [5], will only be summarized 
here. The terms In brackets In Eqs. (35) and (36) can be shown to generate, 
after expansion, terms In sin ^ , cos , sin 3^ and cos 3’l» but the aver- 
aging will retain only the terms In sin if In Eq. (36) and only the terms In 

cos t In Eq. (35), these terms being affected by a coefficient of one-half. 
Thus Eqs. (35) and (36) become: 



(uBG + 


K^0^)/u 


(37) 





(38) 


where 


r(f^w ♦ cg^w) sdg^u / ’’“X '’^^2 * ^ ***®2 “ 

K — j— T — ♦ ~rn *( ~j — rr"5 — + - (^2 ♦ cg^) 

r % r * s \»/ r+s s 


s* 


r(-f,« - cg-w) • sdg.w /rwV r(f, 4 eg.) 4 sdg, w 


T7“T 

r 4 s 


/ r(i)\ r(f^ 4 cg^) 
\s / r^ 4 J 


♦ - (f| 4 cg^) 
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wh«re, from Eqt. (38, 18, 8, 2^) 


f, - -0/I6)«^I3(1 - ♦ b^J + 04)^0 - a)/8 

-Ya^«(U|«/2 - Kb ♦ U2b/2) - youM^ 

-YaU^wd - a)b/2 ♦ (o o^2) («^ + ab^)/8 
+Ya((1 " a)U|W - bwUj] ♦ 0/8)yoiU^^[( 1 - a)^ + (1 - a)b^] 
“2YaU^[(l • a)(K - ^ 2 / 2 ) - bU^/2] 
fj “ ouj^d - a)b/8 - 3a<i)^b/8 - >a^u[a(U2/2 + K) - bUj/2] 

-YOu)(2K + U 2) + YaU^wId - ^)^/k + 3b^/i»] 

+ (1/8) (a + a^/2)(a^b ♦ b^) + Ya((l - a)U 2 tu + bU^uj] 

-(l/8)YaO^ [d - a)^b + b^] - 2 yoU^[( 1 - a)Uj/2 - b(K + U2/21 

g, - I/I6 - ywU^/2 - (l/l|)a«^{a[3d - a)^/k ♦ b^/k] - (1 - a) b^2} 

-d/2)ou)^((I - a)(a^ + 3b^)/4 + ab^/2] 

92 - -yw(K + U 2 / 2 ) + d/^>au^{ab(] - a)/2 - b[(I - a)^/k + }b^/k]} 
(1/2)ouu^[abd - a)/2 + b(3a^ + b^)/<i] 
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Only B In Eqs. (37) and (38) needs to be determined; we can either solve 
numerically for the roots of the fourth-order polynomial given by equating the 
determinant of Eq. (10) to zero or, since B will be a small quantity, we can 
solve analytically the polynomial neglecting terms of order two and higher In 
B . This second method Is presented In [5]. 

STABILITY ANALYSIS 

The effect of noni Inearl ties In t^s critical region can now be Investi- 
gated from Eqs. (37) and (38). Equation (37) Indicates that a positive 6 is 
destabilizing since It leads to > 0 . Also, If , which represents 
the effect of nonlinearities, is positive, we have also a destabilizing effect. 

The only way the nonI Inar I ties can change a stable behavior into an un- 
stable one, and conversely. Is for B and K, to be of opposite sign. When 

this Is the case there Is an amplitude 0, . for which 6 *0 and which Is 

L.C. av 

known as the amplitude of the limit cycle (L.C.). The L.C. Is said to be stable 
If for 0 greater or smaller than 0j^ ^ the amplitude of the motion approaches 
0|^ ^ . The L.C. Is said to be unstable If for 0 greater or smaller than 0j^ ^ 

the amplitude moves away from 0^ 

Hence, from the previous discussion, we can Identify three different pos- 
sibilities, as far as stability is concerned: a) < 0 and 6 > 0 which 

leads to a stable limit cycle; b) Kj > 0 and B < 0 which leads to an un- 
stable limit cycle; c) Kj and B do not meet any of the two previous re- 
quirements. In this last case, 6 and Kj have the same sign and reinforce 
one another, leading to the same conclusion regarding stability as a linear 
analysis. In a) and b) above, the amplitude of the limit cycle is given by: 
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« r “B 

®L.C. Y‘ K, 


( 39 ) 


Figure h shows a plot of ^ as a function of y , the ratio of the 
mass of the fluid to the mass of the pipe plus fluid. Also indicated on Fig. it 
is the type of limit cycle. There are two regions of interest. The first one 
extends from y • 0 to y ■ 1* characterized by destabilizing contri- 

bution of the nonlinearities which thus permit only limit cycles for flow ve- 
locities lower than critical. The second region, from y - .1 to y - .17, is 
characterized by stabilizing nonlinearities and consequent i y limit cycles exist 
only for velocities higher than critical. In the first region, if 0 < 0^^ ^ , 
the oscillation will die out but if 0 > 0^^ ^ , 0 will grow continuously. In 
the second region, if 0 < 0j^ ^ , the motion will increase until it reaches the 
limit cycle and if 0 > 0, > , motion wi 1 1 decay unti I it reaches the limit cycle. 
Nr.turally, the most desirable region is the second one since It permits to operate 
with complete confidence up to th critical velocity and even If one overshoots 
the critical velocity a little, the amplitude of motion will remain bounded. 

For y greater than approximately .17 flutter cannot exist, since ^ 
can no longer be pure imaginary. We also observe that for y close to .1, 
the amplitude of the limit cycle goes to Infinity. This Indicates that the 
nonlinearities contribute nothing to the motion for such a value of y . 

It must be pointed out that in all the above discussions Eq. (38) was 

never mentioned. The reason for this omission lies in the fact that i Is 

av 

the change in the frequency of the motion, which I s of no interest as far as 
stability is concerned. 

EXPERIMENT 

The experimental apparatus (with copper tubes) used by Bohn and Herrmann [3] 
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proved very useful In carrying out tKe following qualitative experiment: 

With the pipe ver'lcally at rest, the fluid velocity Is progressively Increased 
until the slightest disturbance causes the system to flutter. The flow Is 
then reduced by a very small amount and the system Is subjected to small dis- 
turbances. It Is observed that as long as the disturbances are small, the 
motion damps out but that If the Initial conditions are sufficiently large, 
the amplitude of motion grows until It reaches a very large value. This be- 
havior shows the existence of an unstable limit cycle which Is also what the 
theory predicts (Fig. k) since the effective y for the system is about .072. 


CONCLUSIONS 

This study has shown that the effect of finite amplitudes can signifi- 
cantly change the qualitative behavior of the articulated pipe system. Never- 
theless, the results permit us to Increase our confidence in the predictions 
of the linear analysis since Fig. 4 clearly Illustrates that an amplitude of 
motion of at least ten degrees is required to change the flow velocity by only 
a few percent. It thus seems that In most engineering applications the effect 
of the nonllnearl ties can be safely neglected knowing that designers are not 
likely to allow a system to operate within a few percent of the critical 
velocity. 
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